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Abstract 

We study a two-level model coupled to the electromagnetic vacuum and to an 
external classic electric field with fixed frequency. The amplitude of the external 
electric field is supposed to vary very slow in time. Garrison and Wright [Phys. 
Lett. A128 (1988) 177] used the non-hermitian Hamiltonian approach to study 
the adiabatic limit of this model and obtained that the probability of this two- 
level system to be in its upper level has an imaginary geometric phase. Using the 
master equation for describing the time evolution of the two-level system we obtain 
that the imaginary phase due to dissipative effects is time dependent, in opposition 
to Garrison and Wright result. The present results show that the non-hermitian 
hamiltonian method should not be used to discuss the nature of the imaginary 
phases in open systems. 
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1 Introduction 



The natural approach to study quantum systems in contact with an environment is the 
density operator method. It allows one to study quantum states evolving from pure to 
mixed states. In general, exponential decay in the time evolution of the population at 
each excited state of the quantum system occurs due to the exchange between the system 
and its environment. An alternative approach to obtain those exponential decays in the 
probability is through complex energies obtained from a phenomenological non-hermitian 
HamiltonianPl 12 El II]- The non-hermitian Hamiltonian method has been very fruitful 
in approaching various physical problems as, for example, the multiphoton ionization^ 
El in El and free-electron laser theory [HI [7|. As mentioned by Baker in reference jSl, it is 
expected that the non-hermitian Schrodinger equation be a bona fide description of the 
interaction among the parts of a system when the intervals of time are sufficiently short 
such that the coherence state of a subsystem is not destroyed by its interaction with the 
environment. 

Since the disclosure of geometrical phases by M.V. Berry in 19841H1 in cyclic Hamil- 
tonians evolving adiabaticaly, there has been significant search for geometrical phases in 
other physical contexts. For example, Joye et a/. [01 and Berry^O] independently showed 
that the transition probability of instantaneous eigenstates of non-real Hamiltonians in 
the non-adiabatic regime gets an imaginary geometric phase. This imaginary geometric 
phase was measured by Zwanziger et al. in a two-level system [TT]. A great part of the 
work done on geometric phases has been on pure states. Since the work by Uhlmann in 
1986 [121 EHl, however, the study of holonomy has been extended to mixed states under 
unitary evolution[Tll [THl UHl- More recently, Ericsson et a/.[T71 obtained the expression 
of the geometric phase of a quantum system interacting with its environment when the 
unitary evolution of the whole system (including the environment) is known. Certainly 
the discussion of holonomy in mixed states is a very interesting point, since the correct 
expression of the geometric phase for mixed states under unitary evolution is still under 
debate, as in a very recent work by Singh et al.^T^. However, that is not the issue of 
the present communication. Another equally interesting question about the imaginary 
phases in transition probabilities of open systems is the correctness of the application 
of the non-hermitian hamiltonian approach |19t [201 1^ 122] to discuss the nature of such 
phases. 

In particular Garrison and Wright [T^ used the non-hermitian Hamiltonian method to 
study a two-level model with linewidths in the presence of an external electric field, ob- 
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taining the adiabatic limit of the probabihty of this system being in its upper level, after 
the external field has returned to its original configuration. They concluded that the de- 
caying factor has an imaginary correction to Berry's phases. Their result for the two-level 
model coupled to an external classical electromagnetic field with spontaneous emissions 
is already contained in the phases (75) and (76) of reference [3] for any open system de- 
scribed by a non-hermitian Hamiltonian. It is the mathematical structure of these phases 
in the non-hermitian Hamiltonian approach that makes the imaginary phase derived by 
Garrison and Wright in reference ^H] to have an imaginary geometric contribution. 

Garrison and Wright mentioned in their conclusion of reference that the nature 
(time-dependent or path-dependent) of their imaginary phase due to dissipation effects 
should be reexamined using the density matrix approach; this is the aim of the present 
letter. 

In reference [23 studied the adiabatic limit of any periodic non-degenerate Hamil- 
tonian using the density matrix approach (extending the discussion carried out by Born 
and Fock in reference [21] to the density matrix in the basis of instantaneous eigenstates 
of the Hamiltonian), and we concluded that for a quantum system to get an imaginary 
correction to Berry's phases in dissipative phenomena, the functions in the integrals of the 
decay exponentials in the entries of the density matrix would have to satisfy special con- 
ditions. Let e~ -fo ^*'^^^"> be a typical decreasing exponential originated by the presence of 
the dissipative effects. For the integral in the exponential be written as a path-dependent 
integral, the function c{t) must have the form 



and the functions ^E'j(t) have to satisfy two conditions: i) they must not be explicitly 
time-dependent; ii) the time-dependence of functions \I'i(t) must come only from their 



there is no restriction to the regime of the time variation of the set of parameters \<.{t) 
and that it has not to be a periodic function in time, i.e., the path in the k-parameter 
space has not to be closed (see details in reference [2S1)- 

In reference jTHj, Garrison and Wright considered a two-level system interacting with 
a classic external electromagnetic field 





dependence on the set of parameters k(t) = {ki(t), k2(t), - ■ ■ , 



ki(t)). We point out that 



E,ias{t) = Re e£{t)e' 



ii/t 
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where the amphtude S{t) varies very slowly. The two energy states were supposed to have 
a linewidth. They used the time-dependent non-hermit ian Schrodinger (Bethe-Lamb) 
equation in the rotating wave approximation (RWA) ^01 HHj 



dt V C,{t) )-\ Ve-^^' -1^7. ) V C,{t) ) ' ^^^^ 

where 7a and 7?, are the decay rates for the upper and lower levels, respectively, and Caif) 
and Chit) are the corresponding probability amphtude of states \a) and \h). Moreover, 
/S. = uj — u, uj = [Ea — Eb)/h, V = jl ■ e £{t)/2h and Jl = flab is the electric dipole matrix 
element. Letting S{t) = Soe''^^^\ where 0(0) = and 0(T) = 27i, with T > -^^^ they 
obtained a complex Berry's phase /5_ in the expression for the probability of the quantum 
system to be in the state a at time T, 

= -[27r(l - cos(^^o))] = -[27r(l - , ^ ~ )] (2c) 

where 5 = (7a — 76)/2 and Vq = /T ■ e £Q/2h. 

In the present work we want to verify if the result about the path-dependence of the 
imaginary phases derived by non-hermitian method is faithful. To do so, we apply the 
density matrix approach presented in reference jSHI to study the same physical problem 
considered in reference with 7;, = 0. We wish to compare the nature (path or time- 
dependent) of the imaginary phases in the expression of the probability of the atomic 
electron, initially at state |a), being at the same state at time t = T, derived by the 
density matrix approach with the one obtained in the non-hermitian phenomenological 
hamiltonian framework. 

Let the two-level system that describes an atomic electron be represented by the 
Hamiltonian Hg. The atomic electron interacts with a classic external electromagnetic 
field Eciasit) (see eq.(j2i|)) and with the electromagnetic vacuum (the vacuum electromag- 
netic field operator being represented by Eo(a;)). The total Hamiltonian for this model in 
the Schrodinger picture isj^H] 

Ht = He + H/ + Hj„t, (3) 

with 
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k A=l 

and 



He = ^ + Vir), (4a) 
H/ = EE^a;(4a^.^a + ^), (4b) 



Hint = -e(^,i,,(t) + Eo(x)) ■?, (4c) 



where P is the momentum operator associated to the atomic electron, V{r) is the spherical 
interaction potential between the electron and the rest of atom and H j is the Hamiltonian 
of the electromagnetic energy operator of the electromagnetic vacuum. The operator at^ 
(ag_^) creates (destroys) a photon with momentum k in the polarization state A. We also 



have: Uj: = c\k\. 



Following reference p!^, the classic external electromagnetic field Ecias{t) is given by 
eq.()2a|) assuming that S{t) varies very slowly and e is constant. The electromagnetic field 
operators in the vacuum Eo(x) in the Schrodinger picture have the expansion 



k A=i ' 

The dynamics of the density operator of the complete system is driven by the Liouville 
von-Neumann equation. Taking the trace over the electromagnetic degrees of freedom we 
obtain the master equation for the reduced density matrix of the atomic electron p{t) 
written in the basis of the eigenstates of He (He|i) = Ei\i)). The dynamics of p{t) in the 
electric dipole approximation and in the RWA isj2Zll2H] 

at I in 

+7(2a_p(t)a+ - {a^a_,p(t)}) + ^^^l^^f^j^ _ 2A{t)-p{t) 

+2A{t)(j^p{t)a, + Bit) [[ay, p{t)], a,] ) , (5) 



where and a± are the Pauli matrices with a± = ^{a^ ± icxy), pab is the electric dipole 
matrix (we are supposing to be real) with pab = e(a|r|6),(e > 0), 



2 „ 
A=l 



^fc), (6a) 
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being g{k, X) = -ie{k, A) ■ flaty ^e^v 
integration, 



and rj{k) is the density of states introduced in the 



UJQ = 



Ea ~ Eb _ u 

~2h ~ 2' 



(6b) 



^-jr(Ei-Ej)t 



Py(0), 



and 



^ _ 1 



— + 1 



a, 6, 



(6c) 



(6d) 



where ojc is the cutoff frequency that preserves the dipole approximation ( ojc < c/ao, with 
being the atomic Bohr radius). The cutoff Uc becomes a parameter of the effective 
model to be determined by fitting to experimental data. The term corresponds to 
the frequency shift |2ZI- The elements Pij{0) are the initial values of the entries of matrix 
p(t). The functions A{t) and B{t) appearing in eq.(jni) are defined as: A(t) = jldt'piab ■ 
E,Ut') cos[fc^(t - t')] and B{t) = fl^dt'flab ■ E^Ut') sin[^^^(t - t')]. 

We remind that the density matrix of a two-level model must satisfy two conditions: 
i) Tr{p(t)) = 1 and ii) Pbaif) = {Pabif))* , where Pahif) = {(^\p(t)\b)- As a consequence of 
those conditions, the density matrix has only two independent elements; we choose paa{t) 
and Pabit) to be such elements. /^From eq.® the time equations for those two elements 



are 



dt 



Paa{t) 



- 27 - Afiab ■ ^da.(t)Refe*(^"-^'')*G(t))l p,,(t) 



and 



+2pab ■ -Ecias(t)Im 



-i{Ea-Et,)t 



Pab{0) + 2pab ■ Eclasit)Re 



J{Ea-Eb)t 



Git) 



(7a) 



d 

J^Pabit) 



( - 2t{iJo + fi+) - 7 - 2flab ■ ^das(t)G'(t)e'(^"-^^)*)pab(t) 



+2flab ■ Eaas{t)G*{t)e-'^''--''^^'pba{t) + ipab ■ E,ias{t){l - 2paa(0)). (7b) 
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The expression of the function G{t) is: G{t) = dt'jiab ■ Ecias{'t')e~^^^''~^''^^' ■ We see from 
eq. ()7b|l that the real and imaginary parts of pab{t) are coupled due to non-linear effects in 
the classic external electromagnetic field. From eqs.dZ} on, we will be using natural units 

{h = c = i). 

In the same manner as in references j2Sl 12^]; in order to study the T —>■ oo limit of 
eqs.O we apply the transformation 



Pijii), 



i,j = a,b 



and change the time scale to s = t/T; in the limit T oo, one obtains 



d 



T 



27 - 2/i„fc ■ e£o Re{G{s)e' 



-i(<f>{s)-ATs)^ 



(0(s)-ATs)^ 



(9a) 



and 



+flat ■ eSo G*(s)e-^(^W-^^^)rp,,(s) + ^T/i„, ■ e ^(1 - 2p,,(0))e-^(<^(^)-^^^). 

(9b) 

As we write eqs.®, we are assuming that A = 2uJo — ~ The terms proportional to 
^±i{{2a}o+j^)Ts+<f>{s)) contribute in the limit T 00 (as shown in references (221123) 

and the terms proportional to e^''*'^^^'^"''^^'*)) contribute to the dynamics of the density 
matrix in the resonance region when A ~ The function G{s) appearing in eqs.® is 
given by G{s) = flab ■ e ^ /^'^ 

Equation (jUaj) gives us the probability that the atomic electron be at state \a). Its 
solution after one period T is 

PaaiT) = {p..(0) + .^^rft'.(^p,.(0)-^p.,(0))e^(-*'+l^(*')l^) 

+ I ^^.g2,t'rfL_Ml -2(,T+|G(T)P)\ ^IQ) 

2 Jo dt' J J 

Since the function G(t) is an explicit time-dependent function, the condition i in eq.((H) is 
not satisfied and none of the integrals in eq. (|Tn|l is a time- independent integral. We have 
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an overall exponential decay (imaginary factor), but from the definition of the constant 
7 (see eq. (|Haj) ). we obtain that the exponential e~'^"'* is not a geometric (path-dependent) 
imaginary phase, but a time- dependent one as well as the contribution to the decreasing 
exponential coming from |G(T)p. 

For the sake of completeness, we should also examine the solution of eq. ()9b|) after a 
period T. Differently from reference j2H]; here the dynamical equation of Pab{t) couples 
its real and imaginary parts under the regime of strong classic external electric field. We 
get a SU{2) structure for the solution of Pabif)- Calling pi(t) = Pabif) and p2{t) = Pbait), 
the solution of eq. (j9b|l is 



PiiT) 



ij 



dt' 



T e-lT 



J°dt" B(t")-S 



X 



JK 



^g(7t'+Re2(G(t'))~W(G(t')))^^(^^'^ _^ Pj(0)| e(-l)'2ia.oTg-(7T+Re2(G(T))-W(G(T)))^ 



with / = 1,2 and J, = 1,2. In these two last indices we are using the implicit sum 
notation on the r.h.s. of eq. ()ll|) . The symbol T means the time-ordering integrals |3Uj. 
and a are the Pauli matrices. The elements of the column dj{t) are: di{t) = d^it) = ipab ■ 
ef [l-2paa(0)]e-^(<^W-^*). The components of the vector if are: B^{t) = ^^^^P^, By{t) = 



Im(G(t))-^(^-f^»)-Re(6(t))^(^'"f^^^) 



andi3,(t) = -2t [Re(G(t))Im(G(t))] 
We have again the function G{t) that is explicity time-dependent and then the time- 
ordering integrals cannot be converted to path-ordering integrals. Therefore the integrals 
and the phases on the r.h.s. of eq.(fTT|) are all time-dependent, as well as the functions 
with the overall decreasing exponential. 

In reference |E| Garrison and Wright present in eq.(4.9) the value of the geometric 
imaginary phase in the limit of weak electric field and 75 ^ 7a. In the weak electric field 
limit, the SU{2) structure in eq. (fTT|l disappears and the solutions of eqs. lfTIH) and (fTT|) 
become simpler. 



Paa{t) = \ Paa{0) + i 



* f 



,i(At'-</,(t')) 



P6a(0) 



■Ht')) 



Pafe(O) 



-2jt 



(12a) 
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and 



(12b) 

Equations ()12|) have terms with exponential decay (imaginary factors) but they are time- 
dependent, as was obtained in references [2^1122], in opposition to the results of reference 
[in] (see eq.p^). 

Strictly speaking, geometric phases appear in adiabatic processes that happen only in 
the limit T —>■ oo. However, this limit is experimentally implemented by taking T ^ 27i/uj. 
The experimentally measured path-dependent phases are pretty much independent of the 
particular value of T, once the previous inequality is true. In eq.(4.3) of reference [1^] 
we have two imaginary phases: one of them (e"^'^'"^'^-)) does not depend on the value of 
T and consequently is path-dependent, whereas the other one depends on the particular 
value of T. In our results (eqs. ()10|l and (|TT|) ). valid for arbitrary intensities of the classical 
external electric field, all the imaginary phases depend on the chosen value of T. It means 
that the greater the period T of 0(t), the smaller the probability that the electron be in 
the upper level; also, the smaller the correlation pab{t)- 

In summary, we study the nature of the imaginary phase acquired by the probability 
of a two-level model coupled to an external electric field with fixed frequency due to 
its interaction with the electromagnetic vacuum. The interaction among the quantum 
system (two- level system), the environment (vacuum) and the external classical electric 
field is taken in the electric dipole aproximation and in the RWA. Garrison and Wright 
in reference jTHj discussed this same model using the non-hermitian Schrodinger equation 
approach, concluding that the probability has an imaginary phase that is path-dependent. 
In this brief report we apply the matrix density formalism to study its limit of T oo 
and from eqs. ljlUp and (lllj) we conclude that all imaginary phases are time-dependent and 
consequently depend on the chosen value of T. We are considering the situation when 
the effects due to dissipation and time variation of Hamiltonian are of the same order and 
consequently the coherence of initial states are not preserved along the whole period T. 
It is not surprising that those distinct approaches (master equation and non-hermitian 
Hamiltonian) give different results. 

Finally we showed in the present work that even thought the non-hermitian Hamil- 
tonian method has been a very important and useful tool in describing open systems. 
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it should not be applied to discuss the nature (time-dependent or path-dependent) of 
imaginary phases. 
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